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I.  INTRODUCTION 


The  Ballistic  Research  Laboratory,  particularly  the  Advanced  Ballistic 
Concepts  Branch,  has  been  engaged  for  some  tine  In  the  development  of  liquid 
gun  propellants,  one  formulation  of  which  is  a  concentrated  aqueous  solution 
of  hydroxy lamaonium  nitrate  (HAN)  and  various  aliphatic  amine  nitrates.  The 
very  substantial  promise  shown  by  these  propellants  has  recently  been 
discussed*1"5  and  test  fixtures  have  been  successfully  fired.  At  least 
partially  because  pure  HAN  is  difficult  to  prepare  and  handle,  there  has  been 
little  fundamental  research  on  its  thermodynamic  properties.  This  fundamental 
knowledge,  however,  can  reasonably  be  expected  to  become  of  greater  importance 
as  development  of  these  propellants  proceeds.  This  report  will  discuss  the 
state  of  the  current  knowledge  of  the  properties  of  HAN  solutions  and  some 
related  compounds  and  the  possibility  of  computer  simulation  of  them  by 
modeling  them  as  molten  salts.  The  work  that  has  been  completed  on  the 
development  of  computer  codes  for  this  purpose  will  be  summarized,  and 
recommendations  made  concerning  possible  further  work,  both  theoretical  and 
experimental,  on  these  substances. 

II.  THE  MONTE  CARLO  AND  MOLECULAR  DYNAMICS  METHODS 


Simple  monatomic  liquids,  which  can  be  adequately  characterized  by  a 
known  or  assumed  intermolecular  potential  (usually  central),  have  been 
extensively  studied  in  recent  years.  The  most  thoroughly  understood  substance 
is  surely  argon,  and  the  detailed  knowledge  of  the  Intermolecular 
interaction4*5  which  is  available  for  this  simple  substance  is  sufficient  that 
it  appears  that  all  of  its  thermodynamic  properties  in  the  liquid,  gas  and 
solid  phases  can  be  computed  even  to  megabar  pressures,  i.e.  close  to  the 
metallic  transition.  A  comparable  degree  of  knowledge  exists  for  the  other 
rare  gases,  but  there  are  many  substances,  especially  more  complex  ones,  for 
which  the  interactions  are  quite  poorly  characterized. 


When  a  substance  with  a  known  intermolecular  potential  is  studied,  the 
question  of  the  method  of  calculation  of  its  thermodynamic  properties 
arises.  In  a  dilute  gas,  the  low  density  means  that  the  virlal  expansion, 
which  is  a  power  series  in  the  density,  is  rapidly  convergent  and  so  all  the 
equilibrium  thermodynamic  properties  can  be  reasonably  simply  calculated;  the 
transport  properties  of  a  dilute  gas  can  also  be  obtained  from  the  Chapman- 


dW.  F,  Morrison ,  "The  Application  of  Liquid  Propellant  Gun  Technology  to  Field 
Artillery presented  at  the  1982  JANNAF  Propulsion  Meeting, 

^ E .  Freedman ,  ed.,  "An  Outline  of  Liquid  Propellants to  be  published, 

SW,  0 .  Seale ,  "Vulnerability  Studies  for  a  Liquid  Propellant  System ,  " 
presented  at  the  * JANNAF  Propulsion  Systems  Subcommittee  Meeting  (1982), 

*J.  A.  Barker ,  R.  A,  Fisher ,  and  R,  0.  Watts,  Mol,  Phye.,  21,  657  (1971). 

6J.  Barker  and  D .  Henderson,  Revs.  Mod.  Phye,  48 ,  587  (1976). 
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Enskog  and  other  solutions  to  the  Boltzmann  equation.  In  the  other  limit,  a 
solid,  although  its  density  is  high,  has  a  fixed  structure,  and  simple  lattice 
aumsMtions  and  harmonic-  approximation  lattice  dynamics  (and  its 
generalizations)  are  usually  satisfactory.  But  there  are  difficult 
intermediate  cases:  liquids,  dense  gases  and  highly  anharmonlc  solids  (e.g., 
quantum  solids  and  those  near  melting),  where  the  above  simpler  techniques 
fall.  In  these  cases,  although  some  simpler  techniques  such  as  perturbation 
theory  are  sufficiently  accurate  to  be  of  use,  only  machine  simulation 
techniques  will  yield  exact  results. 

The  two  principal  machine  simulation  techniques,  each  of  which  has  some 
variants,  are  the  Monte  Carlo  (MC)  and  molecular  dynamics  (MD)  methods. 

Because  they  have  been  discussed  in  detail  in  many  places,6  only  a  brief 
overview  of  the  usual  procedures  will  be  given  here.  In  each  method,  a  system 
of  a  few  hundred  particles  is  set  up  on  a  computer,  with  the  particles  being 
confined  to  a  cubical  box  whose  side  L  is  chosen  to  reproduce  the  density  of 
the  physical  system  one  wishes  to  simulate.  In  order  to  minimize  surface 
effects,  the  nearest-image  periodic  boundary  condition  convention  Is  used.  In 
the  MC  method  as  formulated  by  Metropolis  et  al.,  particles  make  random 
moves,  which  are  accepted  or  rejected  with  a  probability  proportional  to  the 
Boltzmann  factor  exp(-U/kT),  where  U  is  the  known  or  assumed  potential  energy, 
k  is  the  Boltzmann  constant,  and  T  is  the  Kelvin  temperature.  This  procedure 
is  a  Markov  process,  and  the  average  of  such  equilibrium  thermodynamic 
properties  as  the  energy,  pressure,  heat  capacity,  bulk  modulus  and  radial 
distribution  function  over  the  configurations  generated  with  the  above 
probability  is  a  presumably  exact  calculation  of  these  properties  of  the 
actual  physical  system.  The  Monte  Carlo  method  thus  used  can  provide  all  the 
equilibrium  thermodynamic  properties  of  the  system,  although  it  may  be 
necessary  to  perform  a  number  of  these  simulations  over  a  range  of 
temperatures  and/or  densities  and  Integrate  the  results  with  respect  to 
temperature  or  density  in  order  to  obtain  certain  properties,  especially  the 
entropy. 

In  the  MO  method  due  originally  to  Alder  and  Wainwright,®  the  particles 
move,  not  at  random,  but  in  accordance  with  the  known  or  assumed  force  law 
(the  force  is  the  negative  gradient  of  the  potential  energy  U)*  In  other 
words,  one  simply  Integrates  Newton's  equations  of  motion  for  the  system. 
Rather  than  the  canonical  ensemble  average  used  in  the  Monte  Carlo  method, 
thermodynamic  properties  are  evaluated  as  time  averages  in  the  mlcrocanonlcal 
ensemble,  with  the  agreement  of  the  results  of  the  two  approaches  following 
from  the  ergodlc  theorem  of  statistical  mechanics. 

In  comparing  the  two  techniques,  it  must  be  stated  that  both  yield  all 


«/.  P.  Haneen  and  J.  R.  McDonald ,  Theory  of  Simple  Liquids.  Academia  Press, 
NY  (1976). 

?N.  Metropolis ,  A.  W.  Roeenbluth ,  N.  W.  Roeenbluth,  A.  H.  Teller ,  and  E. 
Teller ,  J.  Chem.  Phye.  21,  1087  (1953). 

J.  Alder  and  T.  E.  Wainuright,  J.  Chem.  Phye.  27,  1208  (1957). 


the  equilibrium  thermodynamic  properties  of  a  system.  Including  the  radial 
distribution  function  (also  called  the  “pair  correlation  function**),  g(r). 

For  purposes  of  comparison  with  experiment,  the  Fourier  transform  of  g(r), 
usually,  called  the  "liquid  structure  factor"  and  denoted  by  S(k),  is  directly 
proportional  to  the  differential  scattering  cross-section  for  elastic  X-ray  or 
neutron  scattering.  The  constant  of  proportionality  is,  however,  different  in 
the  two  cases  because  X-rays  are  scattered  by  atomic  electrons  (Thompson 
scattering),  whereas  neutrons  are  scattered  by  nuclei.  Molecular  dynamics 
yields  also  the  dynamic  structure  factor  S  (k,  <d)  which  is  directly 
proportional  to  the  doubly-differentiated  (with  respect  to  angle  and  energy) 
cross-section  for  inelastic  neutron  scattering.  In  addition,  the  molecular 
dynamics  method  can  yield  some  transport  properties,  especially  the  diffusion 
coefficient,  without  much  difficulty.  Unfortunately  for  the  current 
application,  in  which  the  viscosity  is  of  considerable  interest,  this  is  a 
notoriously  difficult  property  to  calculate  because  it  is  a  many-body  effect 
rather  than  a  single-particle  property.  To  Illustrate  the  difficulty  of  the 
viscosity  calculation,  we  quote  Erpenbeck  and  Wood9  in  discussing  MD  work  on  a 
system  of  hard  spheres  (the  simplest  imaginable  potential)  which  had  consumed 
thirty-four  hours  on  a  Cray  computer:  "...it  is  quite  possible  that  at  high 
density  the  asymptotic  behavior  is  not  established  until  times  substantially 
longer  than  those  attainable  in  the  present  work"  and  "To  obtain  a  better 
estimate  of  (the  viscosity),  then,  would  appear  to  require  both  a  theory  for 
the  decay  of  the  cross  and  potential  terms  and  much  more  extensive  molecular 
dynamics  calculations  for  the  4000-particle  system.  The  latter  would  appear 
to  be  at  the  limit  of  practicality  with  current  computer  speeds." 

The  MC  and  MD  methods  have  been  applied  to  more  complex  liquid  systems, 
e.g.,  those  with  noncentral  intermolecular  interactions  and  systems  with  two 
or  more  different  types  of  molecules.  A  system  composed  of  two  types  of 
molecules  is  characterized  by  three  different  intermolecular  Interactions  and 
three  different  pair-correlated  functions,  the  Fourier  transforms  of  which  are 
called  "partial  structure  factors."  In  the  current  effort,  only  the  MC  method 
has  been  undertaken  to  date,  primarily  because  it  is  easier  to  use.  Until  the 
basic  mechanisms  governing  the  behavior  of  the  liquid  have  been  identified  and 
the  Intermolecular  interactions  well  enough  determined  that  equilibrium 
properties  can  be  reasonably  accurately  calculated,  there  seems  to  be  little 
point  in  attempting  the  more  difficult  calculation  of  transport  properties. 


III.  EXPERIMENTAL  DATA  ON  STRUCTURE  AND  THERMODYNAMIC  PROPERTIES 

The  machine  simulation  techniques  require  as  input  the  Intermolecular 
potentials  and  the  accuracy  of  the  output  calculations  of  thermodynamic 
properties  depends  on  the  accuracy  of  the  potentials.  It  has  become  quite 
plain  Chat  in  the  case  of  HAN  there  are  virtually  no  experimental  data  on  the 
pure  substance  beyond  a  report  by  Ross^  that  its  melting  point  is  48°C. 

There  is  more  Information  available  about  other  hydroxylamnonium  salts.  The 


9J.  J.  Erpenbsek  and  V.  V.  Wood,  J.  Stat.  Phya.  24,  455  (1981). 

1QV.  H .  Robb,  Pros.  Trans.  Nova  Scotia  Inst,  of  Soienee ,11,  95  (1906). 
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structure  of  hydroxy lamine11  itself ,  as  well  as  those  of  hydroxylammonium 
chloride,  hydroxylammonium  bromide,11  and  hydroxylammonium  perchlorate  (HAP) 
are  known.  The  case  of  HAP  illustrates  nicely  the  interplay  between  X-ray  and 
neutron  diffraction.  Dickens,  using  X-ray  techniques  alone,  was  able  to 
make  only  a  partial  structure  determination  because  X-rays,  since  they  are 
scattered  by  atomic  electrons,  do  not  see  protons;  neutrons,  which  are 
scattered  by  nuclei,  can  detect  the  position  of  protons  (or,  more  easily, 
deuterons),  and  the  structure  of  HAP  was  completely  determined  in  a  later 
paper14  by  neutron  diffraction.  The  structures  of  the  solid  phase  of  a 
number  of  nitrate  salts  with  simpler  cations  are  well  known;  for  example,  the 
structures  of  the  nitrates  of  potassium,  rubidium,  silver  and  lithium  were 
reported  in  1913  by  Ewald  and  Hermann.  * 

There  is  no  structure  Information  available  on  the  liquid  phase  of  HAN  or 
on  aqueous  solutions  of  HAN.  There  are  both  theoretical  and  experimental 
determinations  of  the  structure  and  properties  of  some  anhydrous  molten  salts, 
especially  the  alkali  halides.16  There  are  also  neutron  diffraction  structure 
determinations17  of  a  number  of  concentrated  aqueous  solutions,  including 
alkali  halides  and  calcium  and  nickel  chloride.  There  is  no  structure 
information  on  the  pure  liquid  phases  of  other  hydroxylammonium  salts,  but 
neutron  measurements  of  the  structure  factor  of  the  molten  nitrates  of 
lithium,  silver,  sodium,  potassium,  rubidium  and  cesium  have  been  reported  in 
the  literature.16 


IV.  PRESENT  STATE  OF  MONTE  CARLO  WORK 

The  work  done  this  summer  has  provided  a  substantial  start  on  the 
computational  capability  that  will  be  required  for  a  fundamental  understanding 
of  the  equilibrium  properties  of  HAN.  The  specific  accomplishments  are 
threefold.  First,  a  formula  not  previously  found  in  the  literature  has  been 
derived  in  order  to  permit  evaluation  of  the  heat  capacity  at  constant 
pressure  of  a  system  at  a  given  temperature  and  density  by  means  of  a  single 
MC  or  MD  run.  Because  of  its  complexity,  this  formula  will  be  put  into  the  MC 
code  at  a  later  date.  Second,  an  existing  MC  code,  which  due  to  its 


E.  A.  Meyers  and  W.  N.  Lipscomb ,  Acta  Crust,  8,  583  (1955). 

1ZB.  Jerslev,  Acta  Crust.  1,  21  (1947). 

13B.  Dickens ,  Acta  Cryet.  B25 ,  1875  (1969). 

14E.  Prince ,  B.  Dickens,  and  J.  Rush,  Acta  Cryst.  B30,  1167  (1974). 

1SP.  P.  Ewald  and  C.  Hermann,  Z.  Fuer  Krist.  1,  3  (1913). 

ieSee,  for  example,  D.  L.  Price  and  J.  R.  Copley,  Phys.  Rev.  All,  2124  (1978). 

17 J.  E.  Enderby,  in  Ionic  Liquids ,  D.  Inman  and  D.  G.  Loevring,  eds..  Plenum 
Press,  New  York  (1981 1. 

18K.  Suzuki  and  Y.  Fukushima,  Z.  Naturforseh.  32A,  1438  (1977). 
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flexibility  has  been  applied  by  the  author  to  a  wide  variety  of  one-coaponent 
ays teas,  has  been  generalised  to  calculate  the  properties  of  a  two— component 
system  with  arbitrary  short-range  central  Interactions.  Third,  the  MC  code  has 
also  been  generalized  to  study  a  aolten  salt  using  a  suggestion  of  Ceperley 
and  Chester.19  It  will  now  avoid  the  cuabersoae  Ewald  summation  that  is 
usually  required  in  the  case  of  the  Coulomb  potential;  the  code  will  now  do 
(at  least  in  debugging  runs)  only  that  portion  which  must  be  done  by  MC,  the 
remainder  to  be  treated  as  a  perturbation. 


V.  RECOMMENDATIONS  FOR  FURTHER  STUDIES 

Further  progress  on  understanding  the  fundamental  properties  of  these 
systems  will  require  additional  theoretical  calculations  and  experimental 
data,  since  there  has  been  practically  no  work  done  on  them  to  date.  If  it  is 
possible  in  spite  of  the  difficulty  of  working  with  pure  HAN,  it  would  be 
highly  desirable  to  have  structure  information  on  both  pure  anhydrous  HAN  and 
its  aqueous  solutions  as  a  function  of  concentration  from  both  X-ray  and 
neutron  diffraction  experiments.  Any  equation  of  state  and  latent  heat  data 
available  would  be  similarly  helpful. 

Some  knowledge  of  the  three  interionic  Interactions  (those  between  pairs 
of  hydroxylammonlum  ions,  between  nitrate  ions,  and  between  a  nitrate  and  a 
hydroxy laamonium  ion)  is  essential.  The  above  structure  data  will  be  a 
testing  ground  for  these  potentials.  The  most  likely  source  of  initial 
Information  about  these  interlonlc  potentials  is  jib  initio  quantum  chemistry 
calculations.  Although  the  first  step  must  of  necessity  be  to  obtain  and  use 
effective  (spherically  averaged)  central  potentials,  these  potentials  cannot 
be  expected  to  reproduce  the  finer  details  of  these  complex  systems.  As  more 
Information  becomes  available,  a  generalization  of  the  MC  code  to  deal  with 
non-central  interactions  will  be  justified.  In  addition,  electrical 
conductivity,  nuclear  magnetic  resonance,  and  Raman  scattering  measurements 
will  provide  more  knowledge  of  the  basic  mechanisms  that  determine  the 
behavior  of  the  system.  By  shedding  light  on  such  questions  as  the  structure 
of  molecules  within  the  liquids,  solvation,  hydration  and  hydrogen  bonding, 
these  experiments  will  help  to  determine  the  validity  of  modeling  these 
systems  as  molten  salts  and  suggest  how  to  do  more  detailed  (presumably  MD) 
calculations.  As  this  Information  about  the  basic  physical  processes  in  the 
liquids  becomes  available  and  as  the  interactions  are  refined,  the  use  of  the 
MD  code  to  calculate  such  transport  properties  as  diffusivlty  and  electrical 
conductivity  (and  the  viscosity,  provided  there  is  some  advance  in 
computational  techniques)  will  make  sense. 


19D.  M.  Ceperley  and  G .  V.  Chester,  Phye.  Rev.  A1S,  755  (1977) . 

9 


VI.  ACKNOWLEDGEMENTS 


I  would  like  to  express  ay  sincere  thanks  to  the  U.  S.  Army  Research 
Office  for  the  summer  appointment  during  which  this  work  was  done  and  to  the 
Ballistic  Research  Laboratory  and  to  Dr.  Ingo  May  and  the  Advanced  Ballistic 
Concepts  Branch  for  their  hospitality.  Useful  conversations  with  Dr.  Walter 
F.  Morrison,  Jr.,  Dr.  Nathan  Klein,  and  Dr.  George  Adams  and  the  members  of 
his  group  contributed  greatly  to  this  work.  Finally,  it  is  a  pleasure  to 
acknowledge  the  great  debt  I  owe  to  Dr.  Eli  Freedman  for  his  assistance  and 
counsel  on  occasions  too  numerous  to  enumerate. 


REFERENCES 


1.  W.  F.  Morrison,  "The  Application  of  Liquid  Propellant  Gun  Technology  to 
Field  Artillery,"  presented  at  the  1982  JANNAF  Propulsion  Meeting. 

2.  E.  Freedman,  ed.,  "An  Outline  of  Liquid  Propellants,"  to  be  published 

3.  W.  0.  Seals,  "Vulnerability  Studies  for  a  Liquid  Propellant  System," 
presented  at  the  JANNAF  Propulsion  Systems  Subcommittee  Meeting  (1983). 

4.  J.  A.  Barker,  R.  A.  Fisher,  and  R.  0.  Watts,  Mol.  Phys.,  21,  657  (1971) 

5.  J.  A.  Barker  and  D.  Henderson,  Revs.  Mod.  Phys.  48,  587  (1976). 

6.  J.  P.  Hansen  and  J.  R.  McDonald,  Theory  of  Simple  Liquids ,  Academic 
Press,  NY  (1976). 

7.  N.  Metropolis,  A.  W.  Rosenbluth,  N.  W.  Rosenbluth,  A.  H.  Teller,  and  E. 
Teller,  J.  Chem.  Phys.  21,  1087  (1953). 

8.  B.  J.  Alder  and  T.  E.  Wainwright,  J.  Chem.  Phys.  27,  1208  (1957). 

9.  J.  J.  Erpenbeck  and  W.  W.  Wood,  J.  Stat.  Phys.  24,  455  (1981). 

10.  W.  H.  Ross,  Proc.  Trans.  Nova  Scotia  Inst,  of  Science  11,  95  (1906). 

11.  E.  A.  Meyers  and  W.  N.  Lipscomb,  Acta  Cryst.  8,  583  (1955). 

12.  B.  Jerslev,  Acta  Cryst.  1,  21  (1947). 

13.  B.  Dickens,  Acta  Cryst.  B25,  1875  (1969). 

14.  E.  Prince,  B.  Dickens,  and  J.  Rush,  Acta  Cryst.  B30,  1167  (1974). 

15.  P.  P.  Ewald  and  C.  Hermann,  Z.  Fuer  Krist.  1,  3  (1913). 

16.  See,  for  example,  D.  L.  Price  and  J.  R.  Copley,  Phys.  Rev.  All,  2124 

(1978).  —  ' 

17.  J,  E.  Enderby,  in  Ionic  Liquids,  D.  Inman  and  D.  G.  Loevring,  eds.. 
Plenum  Press,  New  York  (1981). 

18.  K.  Suzuhl  and  Y.  Fukushima,  Z.  Naturforsch.  32A,  1438  (1977). 

19.  D.  M.  Ceperley  and  G.  V.  Chester,  Phys.  Rev.  A15,  755  (1977). 


11 


DISTRIBUTION  LIST 


No.  of 

No.  of 

Copies 

Organization 

Copies 

Organization 

12 

Commander 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 

Cameron  Station 

Alexandria,  VA  22314 

1 

Commander 

US  Army  Armament  Mi 
Chemical  Command 
ATTN:  AMSMC  -LEP-L 
Rock  Island,  1L  61 

1 

Director 

Defense  Advanced  Research 
Projects  Agency 

ATTN:  H.  Fair 

1400  Wilson  Boulevard 

Arlington,  VA  22209 

1 

Commander 

US  Army  Aviation  Re 
and  Development  C 
ATTN:  AMSAV-E 

4300  Goodfellow  Bl\ 
St.  Louis,  MO  6312 

1 

HQDA  (DAMA,  C.  Church) 
Washington,  DC  20310 

1 

Director 

US  Army  Air  Mobilit 
and  Development  ( 

1 

Commander 

US  Army  Hateriel 

Command 

Ames  Research  Cente 
Moffett  Field,  CA 

ATTN:  AMCDRA-ST 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1 

Commander 

US  Army  Communicatl 
Rsch  and  Developn 
ATTN:  AMSEL-ATDD 

2 

Commander 

US  Army  Armament 

Fort  Monmouth,  NJ 

R  &  D  Center 

ATTN :  SMCAR  -TSS 

SMCAR-SCA,  M.  Devine 
Dover,  NJ  07801 

1 

Commander 

US  Army  Electronic: 

and  Development  < 
Technical  Support  j 
ATTN:  AMDSD-L 

5 

Commander 

US  Army  Armament 

Fort  Monmouth,  NJ 

R  &  D  Center 

ATTN:  SMCAR-LCA,  D.  Downs 

A.  Beardell 
SMCAR -LCE,  N.  Slagg 

SMCAR -LCS,  W.  Quine 

SMCAR -TDC,  D.  Gyorog 

1 

Commander 

US  Army  Harry  Diami 
ATTN;  DELHD-TA-L 
2800  Powder  Mill  R 
Adelphi,  MD  20783 

Dover,  NJ  07801 

1 

Commander 

US  Army  Missile  Co 

2 

Commander 

US  Army  Armament  Munitions 
and  Chemical  Command 

ATTN:  AMSMI—R 
Redstone  Arsenal, 

Benet  Weapons  Laboratory 

ATTN:  SMCAR  -LCB-TL 

P.  Votis 

Watervliet,  NY  12189 

1 

Commander 

US  Army  Missile  Co 
ATTN:  AMSMI-YDL 
Redstone  Arsenal, 

13 


rsnnnooa.  »  :  KW«Kf 


V  ■  A  .V.V.'.  . 


No.  of 
Copies 


Organization 


DISTRIBUTION  LIST 

No.  of 
Copies 


Commander 

US  Army  Mobility  Equipment 
Rsch  and  Development  Cmd 
ATTN:  AMDME-WC 
AMSME-RZT 

Port  Belvolr,  VA  22060 
Commander 

US  Army  Tank  Automotive 
Command 

ATTN:  AMSTA-TSL 
Warren,  MI  48090 

US  Army  Research  Office  Durham 

ATTN:  R.  Singleton 

P.0.  Box  12211 

Research  Triangle  Park.NC 

27709 

Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL 
White  Sands  Missile  Range 
NM  88002 

Office  of  the  Chief  of 
Naval  Operations 
ATTN:  Code  NOP-351G 
Washington,  DC  20350 

Commander 

Naval  Sea  Systems  Command 
ATTN:  J.W.  Murrln  (SEA-62R2) 
National  Center 
Building  2,  Room  6E08 
Washington,  DC  20360 

Commander 

Naval  Surface  Weapons  Center 
ATTN:  W.C.  Wleland 
Dahlgren,  VA  22448 

Commander 

Naval  SurfaceWeapons  Center 
ATTN:  0.  Dengel 

K.  Thorsted 

Silver  Spring,  MD  20910 


Organization 


Commander 

Naval  Weapons  Center 
ATTN:  C.  Mallory 
S.  Wood 

China  Lake,  CA  93555 
Commander 

Naval  Ordnance  Station 
ATTN:  K.  Mueller 
G.  Poudrler 

Indian  Head,  MD  20640 

Superintendent 

Naval  Postgraduate  School 

Monterey,  CA  93940 

AFATL/ATWG,  0.  Helney 
DLD,  D.  Davis 
Eglin,  AFB,  FL  32542 

AFOSR/NA  (L.  Caveny) 

Bldg.  410 

Bolling  AFB,  DC  20332 
AFWL/SUL 

Kir t land  AFB,  NM  87117 
HQDA 

DAMA-ART-M 
Washington,  DC  20310 

Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905 

Commander 

US  Army  Development  & 
Employment  Agency 
ATTN:  M0DE-TED-SAB 
Fort  Lewis,  WA  98433 

Director 

Jet  Propulsion  Lab 
ATTN :  Tech  Lab 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109 


-.V _V .VW 


No.  of 
Copies 


Organization 


DISTRIBUTION  LIST 

No.  of 
Copies 


Organization 


Director 

National  Aeronautics  and 
Space  Administration 
Manned  Spacecraft  Center 
Houston,  TX  770S8 

The  BDM  Corporation 
ATTN:  Dr.  T.P.  Goddard 
P.0.  Box  2019 
2600  Cearden  Rd 
Monterey,  CA  93940 

Calspan  Corporation 
ATTN:  E.  Fisher 
P.0.  Box  235 
Buffalo,  NY  14221 

Food  &  Machinery  Corporation 
Northern  Ordnance  Division 
ATTN:  J.  Oberg 
Columbia  Heights  Post  Office 
Minneapolis,  MN  55421 

General  Electric  Ordnance  Opt 
ATTN:  J.  Mandzy 
R.E.  Mayer 
H.  West 

100  Plastics  Avenue 
Pittsfidld,  MA  01201 

General  Electric  Company 
Armanent  Systems  Department 
ATTN:  E.  Ashley 
M.  Bulman 

Burlington,  VT  05401 
Director 

US  Bureau  of  Mines 
ATTN:  R.A.  Watson 
4800  Forbes  Street 
Pittsburgh,  PA  15213 

Director 

Los  Alamos  Scientific  Lab 
ATTN:  D.  Butler 
P0  Box  1663 
Los  Alamos,  NM  87544 


Science  Applications,  Inc. 
ATTN:  R.  Edelman 
23146  Cumorah  Crest 
Woodland  Hills,  CA  91364 

AFELM,  The  Rand  Corporation 
ATTN:  Library-D 
1700  Main  Street 
Santa  Monica,  CA  90406 

TRW  Systems 

ATTN:  Rl-1032,  E.  Fishman 

One  Space  Park 

Redondo  Beach,  CA  90278 

Director 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
Johns  Hopkins  Road 
Laurel,  Md  20707 

Director 

Chemical  Propulsion  Info 
Agency 

The  Johns  Hopkins  University 
ATTN:  T.  Christian 
Tech  Lib 

Johns  Hopkins  Road 
Laurel,  MD  20707 

Pennsylvania  State  University 
Dept,  of  Mechnical  Eng 
ATTN:  K.  Kuo 

University  Park,  PA  16802 

Princeton  Combustion  Rsch 
Laboratories,  Inc. 

ATTN:  N.A.  Messina 

M.  Summer field 
475  US  Highway  One 
Monmouth  Junction,  NJ  08852 

Director 

National  Aeronautics  and 
Space  Administration 
ATTN:  MS-603,  Tech  Lib 

MS-86,  Dr.  Povinelli 
21000  Brookpark  Road 
Lewis  Research  Center 
Cleveland,  OH  44135 


tr. 


>»•  v  v*\ 


ft  S' 


DISTRIBUTION  LIST 


No*  of 
Copies 


Organisation 


1  SRI  International 

ATTN:  Code  L3106.  G.A.  Branch 
333  Ravensvood  Avenue 
Menlo  Park,  CA  94023 


1  University  of  Arkansas 
Department  of  Chemical 
Engineering 

227  Engineering  Building 
Fayetteville,  AR  72701 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 
Cdr,  USATECOM 

ATTN:  AMSTE-TO-F 

Cdr,  CRDC,  AMCCOM 
ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-SPS-IL 


16 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comment s/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ _ Date  of  Report _ 

2.  Date  Report  Received  _ _ _ 


3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 
other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


Name 


CURRENT 

ADDRESS 


Organization 


Address 


City,  State,  Zip 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD 

ADDRESS 


Organization 


Address 


City,  State,  Zip 


(Remove  this  sheet  along  the  perforation,  fold  as  indicated,  staple  or  tape 
closed,  and  mail.) 


V  V*  \  \.\  VV  v\«Vv  V  V  V* /  *s  .  V/*  .V  AV-  ;»/'  .SN 

'Sa  a -Vav.v. -V*.  *.  a  V  .  V  \  a  n'w  v’.  n.wv'  v %•  o.w o o 


